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Motile eukaryotic cells propel themselves in viscous fluids by passing waves of bending deformation 
down their flagella. An infinitely long flagellum achieves a hydrodynamically optimal low-Reynolds 
number locomotion when the angle between its local tangent and the swimming direction remains 
constant along its length. Optimal flagella therefore adopt the shape of a helix in three dimensions 
(smooth) and that of a sawtooth in two dimensions (non-smooth). Physically, biological organisms 
(or engineered micro-swimmers) must expend internal energy in order to produce the waves of defor- 
mation responsible for the motion. Here we propose a physically-motivated derivation of the optimal 
flagellum shape. We determine analytically and numerically the shape of the flagellar wave which 
leads to the fastest swimming while minimizing an appropriately-defined energetic expenditure. Our 
novel approach is to define an energy which includes not only the work against the surrounding fiuid, 
but also (1) the energy stored elastically in the bending of the flagellum, (2) the energy stored elas- 
tically in the internal sliding of the polymeric filaments which are responsible for the generation of 
the bending waves (microtubules) , and (3) the viscous dissipation due to the presence of an internal 
fluid. This approach regularizes the optimal sawtooth shape for two-dimensional deformation at 
the expense of a small loss in hydrodynamic efficiency. The optimal waveforms of finite-size fiagella 
are shown to depend upon a competition between rotational motions and bending costs, and we 
observe a surprising bias towards half-integer wave- numbers. Their final hydrodynamic efficiencies 
are above 6%, significantly larger than those of swimming cells, therefore indicating available room 
for further biological tuning. 



I. INTRODUCTION 

The locomotive capabilities of microorganisms are intimately tied to the properties of the surrounding fluid medium 
On scales relevant to most microorganisms, inertial effects are dominated by viscous dissipation; hence, the ejection 
of momentum into the fluid by the shedding of vortices, as observed in the locomotion of fish and birds, is not a viable 
means of propulsion for bacteria and spermatozoa. Instead, microorganisms have evolved to exploit hydrodynamic 
drag. Biological locomotion in this regime is the topic of a vast body of research, and we refer the reader to an excellent 
introduction by Purcell [2] , and the classic texts by Lighthill [3] and Childress 14| . One of the most commonly observed 
means of microorganismic propulsion is the propagation of periodic waves down the length of a slender flagellum. Drag 
anisotropy in viscous flows, in combination with the time-irreversibility of uni-directional beating patterns, renders 
this locomotive form one of rather few relatively efficient means of hydrodynamic propulsion in viscous fluids. 

Due to its ubiquity in Nature, flagellar locomotion has long attracted the attention of biologists, mathematicians, 
and engineers alike. Continuous advances in imaging have revealed new details regarding the structure and kinematics 
of eukaryotic flagella [ES, Q| , but theoretical considerations of flagellar locomotion extend back to the seminal works 
of Taylor 3, 

Hancock |9| , Gray [lO| , and Lighthill [llj . In these works the authors have considered the hydrodynamics 
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of slender body locomotion, developed a resistive force theory for the relationship between velocities and forces, and 
have deduced consequences regarding possible, and in some cases optim al g eometries. Corrections to the simplified 
resistive force theory are found in a more detailed slender body theory 11, 12 .[l3l . [l^ .[l5|. The comparison of theory to 
experiments was furthered significantly in the 1970s in the works of Machin 16|, Higdon \v^ . and Brokaw 0,0, 20 1. 
An excellent review article on flagellar and ciliary propulsion from that era is provided by Brennen and Winet 21 1 . 
More recently, attention has been paid to the relationship between internal structure and hydrodynamics. Camalet 
and Jiilicher 22| have shown that periodic bending and sliding of the microtubule structure of axoncmal flagclla 



can lead to wave generation and organism propulsion. Reidel-Kruse et al. 23| have considered the coordination of 
dynein motors in beating spermatozoa, and have argued that the only theoretical motor coordination that fits their 
experimental data is interdoublet sliding. Other avenues of current active research include the swimming dynamics of 
bodies in non-Newtonian fluid environments, such as the propulsion of spermatozoa in the human female reproductive 
tract 0, 0, 0, 0, H Q ■ 

It is natural to ask about the optimality of the flagellar shapes exhibited by nature. Lighthill Q included a response 
to this question by maximizing a hydrodynamic efficiency over the passage of periodic waves down the length of an 
infinitely long flagellum. He showed that the optimal flagellar shape was one for which the angle between the local 
tangent to the flagellum and the swimming direction was constant. In three dimensions, this leads to an optimal 
flagella in the shape of a rotating helix, a swimming mechanism frequently observed in Nature. In contrast, in two 
dimensions the optimal shape is non-smooth and adopts a sawtooth form. Other early work in this vein was performed 
by Pirroneau & Katz 30| , who considered the hydrodynamically optimal shape of finite slender swimmers, and noted 
an amplitude to wavelength relation for optimality in sawtoothed and small amplitude sinusoidal waveforms. The 
optimal shapes of finite sawtoothed, sinusoidal, and other curves that are amenable to analysis have also been studied 
by Silvester & Holwill [sil, Higdon 17 1, and Dresdner et al. 32 1. More recently, Tam 33] has shown numerically that 
the optimal slender swimmer in a Stokes flow does in fact pass periodic waves along its length, limiting to nearly the 
sawtoothed result of Lighthill. 

In this paper we consider a physically- motivated approach to the question of optimality in planar flagellar locomotion 
by explicitly taking into account the internal nature of the flagellum. We study changes to the hydrodynamically 
optimal but non-smooth shape of Lighthill when internal energetic costs are included. Specifically, we determine 
analytically and numerically the shape of the flagellum which, through the passage of a wave down its length, swims 
the fastest while minimizing a newly-deflned swimming energy. This energy not only includes dissipation in the 
surrounding fluid, but also (1) elastic energy stored in the bending of the flagellum, (2) the Hookean energy stored 
in the relative sliding of the polymeric filaments (microtubules) which create the waves of deformation, and (3) the 
internal dissipation due to the presence of a fluid inside the axoneme. In the case of infinite fiagellum length, we 
show that this approach regularizes the non-smooth solution of Lighthill, at the expense of a small decrease in the 
hydrodynamic efficiency (from 8.5% to 7.5%). Finite- length swimmers also display smooth flagellar shapes, and 
we show that optimal shape is determined by a competition between rotational motions and bending costs with a 
surprising bias towards half-integer wave-numbers. The hydrodynamic efficiencies of the optimal finite-length flagella 
are above 6%, significantly larger than those of biological cells (typically in the 1% range), indicating available room 
for further biological morphological tuning. 

The paper is structured as follows. In f|TI]we introduce the notation for the swimming kinematics, as well as the 
new swimming energy measures we use in the paper. The case of an infinite-length flagellum is treated analytically in 
§1111 bv a variational approach, and solved numerically for finite bending and sliding costs. The results for finite-size 
flagella are presented in piVI We conclude with a discussion of our results and their implications for the biophysics 
of motility in f|Vl 
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II. KINEMATICS, FLUID-BODY INTERACTION, AND ENERGETIC COSTS 

A. Kinematics 

We consider the passage of a periodic waveform down along an inextensible flagellum of length L and radius a, 
which is confined to motion in the x ~ z plane. The waveform is described by X(s) = {X{s), Z{s)), where s G [0, L] 
is the arc-length (|Xs(s)| = 1). The waveform is chosen so that the body is initially oriented along the x-axis, with 
X(0) = 0, and 

X{s + A)=X{s) + X, Z{s + A) = Z{s), (1) 

where A is the distance along the flagellum between wavelengths. A is the physical wavelength, so that a — X/ A < 1 
is a contraction factor due to the waviness of the flagellum [3| . We define L = kA, with k the number of wavelengths 
along the body (not necessarily integral). 

The body motion is illustrated in Fig. [TJ At time t, the waveform {X{s), Z{s)) is assumed to pass along the length 
of the body at an angle 6{t) to the x-axis. In a frame of reference moving with the traveling wave, the flagellum 
moves tangentially with uniform speed c and a period T* — A/c. Defining Xo(t) as the position of the head and 
r(s,t) = X(s — ct) — X(— ct), the body centerline is written as 

x(s,<) =xo(t)+Rr(s,t), (2) 

where i — cos{6)x + sin(^^)z, and R is the rotation operator, 

cos^W - sin 0(0 \ 
~ y s\n9{t) cos 9{t) j ' 

Hence, the velocity of each point may be written in the lab frame as 

u{s,t) = ±o{t)+R(rt + eit)r^) , (4) 



with X^ = {—Z.X) and rj — — cXs(s — ct) + cX.s{—ct). The unit tangent vector along the body (in the direction 
of increasing s) is denoted by s = R(X^), and the normal vector is n = s^. 



B. Fluid-body interactions 



The fluid-body interactions are modeled using the local resistive force theory of Gray & Hancock (1955) [34 1. 
Resistive force theory relates the local fluid force per unit length and the local fluid/body velocity (equivalent by the 
assumption of a no-slip boundary condition). In the classical theory, the velocity of the slender body at a station s is 
separated into tangential and normal components, as is the corresponding force per unit length f (s, t): 

s • f (s, t) — Kt s • u, n • f (s, t) = Km n • u. (5) 

The force per unit length on the fluid, f(s,t), is thus taken to be 



f (s, t) = Ktss^m + Kn{I - ss^)u = {Kt ~ Kn)s{s ■ u) + Knu. (6) 
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FIG. 1: Swimming fiagellum with notation. A periodic waveform with physical wavelength A is passed from the head xo(f) 
down along the body centerline (to the right) at an angle 9{t) to the horizontal. The body moves opposite the direction of the 
wave in mean. Unit tangent and normal vectors are also indicated. 



With no external forcing, the dynamics are thus set by conditions ensuring zero net force and zero net torque on the 
body at all times 

f{s,t)ds = Q, / (x(s,i) -Xo(t)) X f(s,t)ds = 0. (7) 
Jo 

These three equations are linear in the velocities Xq and 6, which are solved by a simple matrix inversion. The 
time-dependent body geometry determines uniquely the velocities at all times. 

Corrections to the resistance coefficients, and to the local theory in general, are the subject of a number of studies 



[ill. l34l. l35l| . For this study wc fix the ratio rj. = Kj^f/Kx = 1/2 with the acknowledgement that this ratio has been 
found in these other works to be dependent (though logarithmically) upon the ratio a/ A which we take to be very 
small. Non-local effects are potentially significant for the study of all but the thinnest of bodies; however, the use 
of slender body theory is complicated by the need for high resolution of the body shape near any regions of rapid 
geometric variation such as a kink. The slender body theory generates a system of Fredholm integral equations of 
the first kind which are in general susceptible to oscillatory behavior or slow convergence in their numerical solution, 



particularly when the immersed boundary has a sharp geometry [36l |. This said, we note that Tam [33.] shows the 



near recovery of Lighthill's sawtoothed waveform solution using the full slender body theory. 

C. Energetics 

1. Dissipation 

The rate of mechanical work done by the body against the fluid, $*(t), is determined through an integration along 
the body centerline, 

^*{t)= [ i{s,t)-u{s,t)ds {>0), (8) 



and is seen to be non-negative due to the form of Eqn. [H] Averaging over one cycle, we define 



$* = (|.*(t)). 



(9) 
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where 

(^*W> = ^^'^ ^*it)dt. (10) 

In addition to performing work against the fluid, internal forces must also be exerted in order to create bending 
waves along the flagellum. In this paper we are considering these forces by explicitly taking into account the elastic 
nature of the flagellum, as well as the presence of internal dissipation. Three new measures of energy are therefore 
defined below. 



2. Bending 



Figure [2] shows a TEM image and a cross-sectional diagram of a typical Eukaryotic flagellum, in this case that of 
the organism Chlamydomonas. The internal structure of a Eukaryotic flagellum, known as the axoneme, is usually 
composed of nine microtubule doublets which encircle a central microtubule pair (though other numbers and mod- 
ifications of this basic pattern have been observed) 2l| . Dynein molecular motors act to generate shear forces that 
cause sliding between the outer doublet microtubules, and consequently the macroscopic passage of waves along the 
flagellar length 37|. Nexin proteins are elastic links that act to keep the outer microtubule doublets well spaced. 




Microtubule doublet 
Nexin 



Inner dynein arm 
Outer dynein arm 



Membrane 



Radial Spokes Central singlet microtubules 



FIG. 2: The structure of the flagellar axoneme is seen in a (a) TEM image (Rippel Electron Microscope Facility, Dartmouth 
College) and (b) a cross-sectional diagram of a Chlamydomonas flagellum. 



We model the elastic energy stored in the bending of the axoneme, SBendingW^ ^ function of the flagellum's 
effective Young's Modulus E, its second moment of inertia /, and the local flagellum curvature K{s,t), 



^BenMngit) = EI K^{s,t)ds, 



(11) 



(see 38|). Kink instabilities have been shown to form when soft elastic cylinders are bent beyond a critical radius of 



curvature 



39|; however, we take such defects to be negligible given the assumption of the vanishingly small aspect 



ratios considered here. A bending power associated with Eqn. (fTTj) is defined as 



'^Bending ~ Bending i^)) ' 



(12) 



and represents the time-averaged elastic energy stored in the flagellum per unit period of the wave. 
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3. Elastic sliding 



The relative sliding between the microtubule doublets (Fig. ^ is understood to account for the generation of 



bending moments and large scale undulations [2^, |23|, |40[. In the study of planar waves, a common abstraction of 



the internal sliding is to consider a "two-dimensional axoneme" as illustrated in Fig. [3] (following Camalet & Jiilicher 
22l|). A rigorous connection between the two-dimensional consideration above and the full three-dimensional axoneme 
is presented in 41], and bending moment propagation in flagella by such sliding action is considered in 42 1. 



T 
a 




X{s) 



FIG. 3: Two-dimensional analogue of the axoneme shown in Fig. [2] A local sliding displacement A is defined as the difference 
between the arc-lengths of the top and bottom curves X(s) ± (a/2)n(s). Here we have set zero relative displacement at the 
leftmost edge (A(s) = 0). 



To capture an energetic cost due to a material shear of this nature we define a sliding energy per wavelength under 
the assumption of a Hookean internal response. Based on the two-dimensional structure described above, we define a 
local sliding displacement A(s,t) as a difference between the arc- lengths of the top and bottom curves (as illustrated 
in Fig. El), 



A(s,i) = A(s-,i)+^' (|a.(x+^n) |-|9,(x-^n) |) 



ds' 



A{s,t)+a J K{s',t)ds' 



(13) 

(14) 
(15) 



Here we have defined the local tangent angle tlj{s,t), where Xg — cos{ip) and Zg = sin('i/;), and tpg = k. A complete 
description of the sliding distribution requires initial specification of the relative distance between top and bottom 
curves at a single reference point A{s = s,t — 0). Subsequently for t > the sliding dynamics are set by the time- 
dependent geometry as we will show. The relative sliding displacement changes with the local shears generated by 
bending moments, and the behavior at s = s may be written in terms of a jump in upper and lower curve velocities, 

At{s,t) = [s • u] as • fit = -(ac) s • = -(ac) ipsisjt), (16) 

where [/] = /*°p — jbottom^ Hence, the sliding displacement may be written as 

A{s,t)^A{s,0)-{ac) f i;s{s,T)dT + a{ij{s,t)-ij{s,t)). (17) 
Jo 
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We exploit the travelling wave structure, ips = ~\4't: and integrate the above to give 

A(s, t) = (A(s, 0) - a V'(s, 0)) + a ^(s, t) (18) 
= a(cA+V'(s,t))- (19) 

Reidel-Kruse et al. [^I set s = and couple the sliding displacement at the flagellar base, A(s — 0,i), to the 
internal sliding dynamics along the body length. A sliding energy is defined through a shear modulus G as 

^*suMna{t)=G f A'ds, (20) 
Jo 

and we define an associated sliding power as the amount of sliding energy stored per unit period of the wave, 

'^Sliding = {^Slidingi^)) 1"^* ■ (21) 

For periodic waveforms, (cca) is equivalent to the period-averaged sliding displacement, and we refer to ca as a 
dimensionless base sliding. Having assumed a Hookean elastic response to sliding, ca 7^ therefore corresponds to 
a non-zero net internal moment. This moment would act, absent any other internal forces, to drive the flagellum 
towards the ca = state, or 



A(s,0) = aV(s,0), (22) 

with the sliding displacement precisely equivalent to the tangent angle multiplied by the body thickness. This special 
case is illustrated in Fig. |4l along with an illustration of a body with non-zero base sliding, ca > 0. For a given body 
shape, the second corresponds to a state with a larger internal energy. Absent external forces there can would also 
be a net body rotation in an energy minimizing response to ca 7^ 0. 




FIG. 4: Two arrangements of sliding displacement density are illustrated. The first corresponds to zero mean sliding displace- 
ment, (aca) = 0. The second corresponds to a positive mean sliding displacement, (oca) > 0, and hence a non-zero net 
internal moment. 



8 



4- Rate of sliding 

The locomotive properties of many organisms may also depend upon the dynamics of a fluid internal to the body. 
Hence, we consider a third cost of locomotion, that of internal dissipation due to the sliding of an upper and lower 
boundary as described above, 

^*Internal (0 ^ ^ = ^A^^ ^ (^3) 

where [ii is the viscosity of the internal fluid and we have used Eqn. 1161 As a result, the optimal waveforms determined 
by including an elastic bending cost will also capture the effects on the optimal shape determined in the presence 
of an internal dissipation cost. This is at first glance a puzzling result; a shear of two surfaces certainly need not 
require a curvature. However, given the assumption of a periodic waveform and the associated velocities, the shears 
generated are in fact due only to the bending of the fllament. Without loss of generality, we therefore consider below 
only elastic bending and sliding costs. 



5. Efficiency 

The system is made dimensionless by scaling velocities on the speed c, lengths on L (A) for flnite-length (infinite- 
length) bodies, and time on the ratio of length to velocity scale. Geometrical variables, velocities, and forces are 
henceforth understood to be dimensionless, and the period of motion is denoted by T. The rate of mechanical work, 
and bending and sliding powers for flnite bodies are written as 



^={j f-uds^, V Bending = 1b{^ j ds^ , Vsiiding = Is J A^ds^, (24) 

where 7^ and 7^ are dimensionless and may be inferred from the deflnitions above. Finally, we define a swimming 
efficiency. A common measure of hydrodynamic efficiency, rjH, for low Reynolds number swimming is the ratio of 
the rate of work required to drag the straightened ffagellum through the fluid to the rate of work done to propel the 
undulating body at the same velocity, 

m = — (25) 

where U = |(xo(t))| is the (dimensionless) mean swimming speed. 

In this paper we define a generalized swimming efficiency, rj, by including as well the rates of energetic expenditure 
due to internal bending (or internal dissipation) and sliding, 

„ = (26) 

il-AB){l-As)^+iAB/lB)VBend^ng + iAs/-/s)'PsHd^ng' 

Ab G [0, 1] and Ag G [0, 1] are dimensionless numbers which allow for variation of the relative importance of bending 
and sliding energetic costs. The optimal swimmer is henceforth defined as that waveform X(s) which maximizes 
the efficiency 77. We note that an alternative measure of efficiency has been defined in Ref. however, in order to 
compare our work most directly to the classical result of Lighthill and other recent works we have chosen the measure 
above. 
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III. BODIES OF INFINITE-LENGTH 



A. A classical result by a variational approach 

We begin by showing that a variational approach yields the classic result due to Lighthill When the body is 
infinitely long (L ^ oo), it is useful to decompose the body velocity into the tangential motion along the waveform 
and a swimming velocity, U = ([/, 0) (as in [^) which we achieve by defining 

U = xo - (ai - s\s=o). (27) 



Since {ax — s|s=o) = for periodic X(s), this wave-frame velocity is equivalent in mean to the head velocity, 
([/) = (io) — U. The velocity of a point on the body in this special case may then be written as u = ([/ + a)x — s, 
and Eqns. [7] reduce to 

1 - (1 - rfc)/3 

a= f X,ds, 13^ f {Xs)^ds, (29) 
zo = -z- s\s^o, 0{t) = 0. (30) 



The rate of mechanical work is given by 

$ = ((7 + af{l - (1 - rk)l3) + rfe(l - 2alfi + a)). (31) 



Defining the slope function g{s) = Zs{s), and setting the variational derivative of the efficiency to zero {Srj = 0) the 
following algebraic equation is generated for g{s): 

( ^2 1 (l-/3)^(l-(l-rfc)/3)4 

Since the absolute slope \g{s) \ is constant we may compute simply the constants P — 1 — = a^. These relationships 
then yield P — 1/(1 + y/rk), and 

g{s) = ±J-^^- (33) 

This is the result of Lighthill For rk — 1/2, the physical slope of the sawtoothed waveform is sin~^(|(7(s)|) = 40.06°. 
The associated efficiency isr] = {l — ^/rk)^ = 0.0858, and the swimming velocity isU = — {I — y/Fk) / {1 + y/Fk) = —.224, 
opposite the direction of the travelling wave. Generally speaking, a body shape which alternates between the positive 
and negative slopes ±|g(s)| at an arbitrary number of points yields the same swimming velocity and efficiency. The 
structure even admits pathologically discontinuous g{s), up to the point at which the resistive force theory becomes 
invalid. This complication is removed in the case of infinite-length by simply requiring that a fundamental periodic 
mode is exhibited over a unit wavelength. 

If the bending cost is included {Ab > 0) but the sliding cost is ignored {As — 0), after some algebra we find that 
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the variational problem leads to the following integro-differential equation for the slope g{s): 

. ( , \ ({'^- AB)rk + {AB/ls)VBendzng\ r, J 

+ " \ (l-(l-r.)mi-/3) j ~ ^ ^ 



+ 2.,(AW7B)(^ „_,;— ^,^, J,(l-.-) = 0. (34) 

The ^ limit is readily seen to be a singular one, which is expected due to the lack of regularity in the sawtoothed 
solution. A numerical optimization for the case of general As > is presented in a later section. However, rewriting 
Eqn. [34]for > with its highest derivative alone on the left hand side makes clear that further differentiations 
introduce no irregularities for \g{s)\ < 1, so that such solutions have g € C°°. 



B. Numerical optimization 

Having already assumed the fluid is modeled by the Stokes equations and since length scales out of the efficiency 
measure completely for a body of infinite length, there is a scale invariance in the determination of the optimal 
waveform. We account for this invariance in the infinite-length case by assuming that the waveform cannot be 
decomposed into smaller periodic forms. In other words, we scale the optimization so that the fundamental shape 
is expressed exactly once in the spatial period s e [0, 1]. This is achieved by considering the following basis for the 
tangent angle i)(^s): 

V'(s) = tan"^ ( ^4~T ) 51 cos(27rns), (35) 

subject to the constraint that Z(s + 1) = Zis). We make a simplifying assumption that -0 € [— tt, tt]. The solutions 
found without enforcing this constraint are consistent with this assumption, but the numerical search procedure can 
become unstable or slow in some cases without its application. The phase of the waveform is irrelevant for bodies of 
infinite-length. Given the tangent angle at each station s, we recover the flagellar shape by integration, 

g(s) = Z,(.s) =sin(7^(s)), X, = cos(V'(s)). (36) 

Numerical optimization is performed using a SQP, Quasi-Newton, line-search method in MATLAB's optimization 
toolbox. The tangent angle is written as a finite sum over the first n* Fourier modes, and the corresponding coefficients 
a„ are determined so that the efficiency r\ is maximized. The optimization routine runs until the line search detects 
a local solution gradient with a relative error tolerance of 10^^^. 

The swimming velocity is determined by discretizing the slender filament in the arc-length s by M uniformly 
distributed points and solving the linear system, Eqn. [T] Since the body is an infinitely long traveling wave, the 
velocity and rate of mechanical work are constant in time, and hence are determined numerically at t = alone. 
Hence, each iteration of the optimization routine requires the creation and inversion of only one algebraic equation 
for i;o(i), since ZQ(t) — 9{t) = in the infinite-length case. Generally, the values {M — 4000, n* = 80) are sufficient 
so that further resolution has a negligible effect on the solution. 

The constants (a, /3) as well as the rates of work ($, VBending) are determined by quadrature in s. Given that their 
integrated arguments are periodic on s G [0, 1], a simple trapezoidal rule yields spectral accuracy. The optimization 
routine was seeded with a variety of initial flagellar shapes to increase the probability that a global maximum of 
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efficiency was achieved. However, given enough spatial resolution the solutions found for the infinite-length case did 
not vary, regardless of the initial guess. In addition, the solutions so found have been verified by insertion into Ean. l34l 

C. Finite bending costs: numerical results 

We now present the optimal shapes of infinite- length bodies with the inclusion of the bending cost. Figure [5] shows 
the optimal waveforms for a sequence of bending parameters, Ab, with no sliding cost {As = 0). In order to best 
compare the shapes, the optimal waveforms are rescaled to the same physical wavelength for presentation. 
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FIG. 5: (color online) Optimal waveforms for an infinite flagella with various bending costs {Ab 7^ 0), in the case of no sliding 
cost {As = 0). For a small bending cost {Ab <^ 1) the optimal waveform is very nearly the analytically derived sawtooth 
function. As the bending cost increases, the shape undergoes its most dramatic change for Ab ~ 10"^, and settles to very 
nearly a sinusoid for Ab ~ 1- 

For Ab — 10^^ the optimal waveform is very nearly the analytically derived sawtooth function. As the bending 
cost increases the shape undergoes its most dramatic change near Ab ^ 10-^, and settles to very nearly a sinusoid for 
Ab = 1; specifically, to the Z{s) w .1159 sin(27rs) - .0017sin(67rs), or Z{x) w 0.1208 sin(27ra;/Q;) + 0.0033 sin(67r2:/a), 
with a ~ 0.85096. The inclusion of bending costs as we propose in this paper therefore effectively regularizes the 
non-smoothness of Lighthill's solution. 

In Figs. [6] and [71 we provide a closer inspection of the optimal shape and its properties for the same range of bending 
parameters Ab- First, and as expected, with the additional cost of bending (increasing Ab), the curvature at the 

apex decreases from its infinite value in the Lighthill solution. The wavelength-normalized curvature at the apex is 

1/2 

further plotted in Fig.[7K. It is not surprising to recover an apex curvature k gs ^ Ag due to the form of Eqn. [34l 
For Ab ^ I, Eqn. [34] is dominated by the algebraic expression of Eqn. [32| outside a boundary layer region where 
Ab 5ss = 0(1). The terms in this equation are of like order in a region of size s ^ \J Ab around the discontinuities in 
the Lighthill solution. 

The third column of Fig. [6] displays the slope at the midpoint of the body, which we find is not monotonic in the 
bending parameter Ab- This angle (absolute value) is further plotted in Fig. [7)d. The absolute slope increases from 
40.06° to approximately 49° at Ab ~ IQ-^'^, then decreases to nearly 44° when the bending costs are prohibitively 
expensive {Ab 1). 

The total and hydrodynamic efficiencies are displayed in Fig. [7]: as a function of the bending parameter Ab - For 
Ab > IQ-^, the efficiency decreases approximately like r\ ^ A^ , which may be predicted given an inspection of the 
efficiency measure 77. The limiting value of the total efficiency is 0.00130. Given the relatively small global change in 
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FIG. 6: (color online) A closer look at the optimal waveform with the inclusion of bending [Ab 0) but no sliding costs 
[As ~ 0), and the departure from the Lighthill result (dashed lines). 



the optimal shape, the hydrodynamic efficiency does not decrease as drastically as the total efficiency with increasing 
bending costs. For Ab — > 1, the hydrodynamic efficiency approaches the limit rjn — 0.0746, a decay of only 13% from 
the optimal Lighthill solution in the case where bending is without cost. In this case, the design of an organism or 
manmade swimmer is far more sensitive to the energetic costs due to bending than to hydrodynamic costs. 

The swimming speed, U, is shown in Fig. [TJi as a function of the bending parameter As, which decreases from the 
analytical solution in the previous section with increasing bending costs. Even though the hydrodynamic efficiency 
only decays 13% from the optimal sawtoothed solution, the swimming speed decreases by approximately 20%. Fig. [7^ 
shows the wavelength-normalized amplitude, bk, where b — 27r||Z||oo and k — 1/a is the wave-number. This amplitude 
decreases from approximately 1.31 to 0.87 as Ab 1- 

Finally, Fig. [7]: shows the odd-numbered Fourier modes in of the optimal waveform for four different bending 
parameters Ab on semi-log axes. We observe that the Fourier coefficients decay linearly on this scale, so that 
a„ ~ iQ-pi^B)n fpj. ig^j.gg (odd-numbered) n and constants p{Ab) as indicated in the figure. This decay, faster 
than polynomial in n, indicates that the optimal solution has g{s) G C°°, in agreement with the previous comment 
(even- numbered Fourier modes are zero to working precision). 



D. Sinusoidal waveforms 

The variational result given by Eqn.[34]does not lend itself to a straightforward long wavelength (small amplitude) 
asymptotic analysis; indeed, the solutions of interest require g{s) = 0(1) so that g{s) ~ g{s)'^ ■ However, a linearization 
of Eqn. [34] balances the elasticity with the hydrodynamics: roughly, gss + uig — 0. Hence, we expect periodic, near- 
sinusoidal solutions in general. This intuition is already corroborated by the near-sinusoidal solution for large bending 
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loglo(^B) 

FIG. 7: (color online) Properties of the optimal waveform for an infinite fiagellum as the bending parameter Ab is varied from 
to 1, with no sliding cost {As ~ 0). (a) Normalized maximum curvature, (b) Absolute midpoint slope (in degrees); the behavior 
here is illustrated in the third column of Fig. |6] (c) Swimming efficiency rj and hydrodynamic efficiency rjH- (d) Swimming 
velocity U. (e) Normalized waveform amplitude bk, with b = 27r||Z||oo and k — 1/a the wave-number, (f) Odd-numbered 
Fourier coefficients |a„| are shown to decay rapidly on a semi-logarithmic scale, indicating the solutions to be infinitely smooth. 
(Even-numbered modes are zero to working precision) . 

parameters {Ab ~ 1), as previously discussed. For comparison, we compute the optimal amplitude of a sinusoidal 
waveform. Inserting the ansatz Z{s) — (6/27r) sin(27rs) into the efficiency measure, Eqn.[2Sl and using Eqns. (|28ll30p . 
yields 

a^-E{b^), P^l-^, VBendrng^^TT^l- Vl-b^), (37) 

where E{m) is the complete elliptic integral of the second kind, (m G [0, 1]). The slope amplitude b that maximizes 
the efficiency is determined by setting db rj — 0. The expression for db t] is unwieldy so a rootfinding algorithm is used 
to locate the optimal amplitude as a function of the bending parameter and the resistance coefficient ratio r^. The 
results are reported in Fig. [H The optimal sinusoid when there are no bending costs {Ab = 0) has an associated 
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efHciency of 0.0782, which is smaher than that of the optimal sawtoothed shape by only 9%. The total efhciency is 
seen to decrease with increasing bending costs to a limiting value of 0.00128 when Ab = 1, which is just barely smaller 
than the efficiency of the fully optimal solution determined by the numerical optimization (rj = 0.00130). Here again, 
the hydrodynamic efficiency does not decrease significantly as the bending becomes more expensive, since the optimal 
shape does not change dramatically. For As — > 1, we find rjH ~ 0.0728. 
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FIG. 8: (color online) Optimal sinusoidal waveforms of an infinite flagellum. (a) The efficiency rj and hydrodynamic efliciency 
r;^ of the optimal sinusoid are shown as functions of the bending parameter Ab- (b) Normafized amplitude bk. (c) Midpoint 
slope. 



The optimal amplitude decreases from a value of 6fc = 1.06 for Ab = to a limiting value of bk — 0.837 when 
Ab = 1. This matches very nearly the previous result in the numerical study {bk = 0.87). Fig.[Hj3 shows the midpoint 
slope angle, which is monotonic in the bending parameter Ab, and decreases from 57° to 45° as Ab — > 1, the latter 
result again nearly matching the result for the fully optimal shape. 



E. Finite internal sliding costs 

We now turn our attention to the consequences of finite internal sliding costs. Given the waveform periodicity 
the time-average of the sliding power may be simplified to an integration against the initial waveform. Recalling the 
travelling wave structure, and using that the tangent angle has zero mean, we have the expression 

VsUMng = 7s + ^ ^ t) ds d^j . (38) 
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The inner integration may be written as independent of time by a simple manipulation, 



1 /.I i<i i<i j-i j-i 

2/'„ +\ j„ jj- _ / / „;,2/'o_/'„ j„ j+ / / „;,2/ 



■\\}\s,t)dsdt^ / / tp-" {2n {s ~ t)) ds dt ^ / / tp'' {2t: s) ds dt (39) 
^0 Jo Jo Jo Jo 



1 fl 



V'n27rs)ds= / V (s,0)ds. (40) 
In other words, the sliding displacement travels with the waveform. In terms of the slope function g{s), 

VsUd^ng = 7s (^4 + ^ [sin-^ (^(s))]^ ds^ . (41) 

We are now prepared to consider the variational derivative of the efficiency when the sliding cost is included {As > 0) 
but the bending cost is ignored {Ab = 0). The result is once again an algebraic relation for the optimal shape: 



ci9{s) + C2g{s)^l-g^s)+C3 sin- \g{s)) = 0, (42) 

where we have defined 

Cl =(1 - /3)(1 - (1 - rfe)/3)2((l - AsH + {As/"^s)VsHd^ng). (43) 

C2 =rfca{(l - As)rk ((a^ - 2) (1 - (1 - r^)/?) + rm'') (44) 

- 2(As/7j(l - (1 - rk)P)Vsl^d^ng]. (45) 

c^=Asa{l~ (3){l^[l~Tk)(3f. (46) 

The form of Eqn. indicates that \g{s)\ = go (constant), so that once again we recover a solution which is not 
smooth at points where the slope switches sign. The inclusion of internal sliding costs does not, therefore, regularize 
the optimal swimming shape. For a sawtooth function we must have P — — \/l — g^- Inserting these values, the 
resulting expression is not amenable to analytical solution, but we apply a simple root-finding algorithm to determine 
the slope go that satisfies Eqn. [42l 

The results are reported in Figs. [H where we show the resultant slope angle -0 = sin~'^(|(7o|) as well as the total 
and hydrodynamic efficiencies for a selection of constants ca- As As 0, the Lighthill solution is recovered in all 
cases: ip ~ 40° and ry = 0.0858. As the sliding cost becomes more dominant we find that the total efficiency drops 
significantly, and the optimal shape is a sawtooth function with a smaller and smaller amplitude. The slope angle 
■0 decreases to nearly 33° for ca = when the sliding costs vastly outweigh the hydrodynamic costs {As 1). 
For non-zero base sliding (ca > 0) the sliding cost is relatively less dependent upon the waveform, and the optimal 
solution does not change as significantly with increasing As- 

The total efficiency decays significantly with increasing As from the Lighthill value r] = 0.0858 for ca = 0, but 
the hydrodynamic efficiency only decreases to a limit 0.0803. Since the optimal shape for very large sliding costs is 
similar to that when there are no sliding costs, it is not surprising to find that the hydrodynamic efficiency only begins 
to noticeably decay for As > 10^^, and even then only to a limit 0.0803. If the base sliding is larger, however, the 
waveform retains a nearly 40° slope, and the hydrodynamic efficiency stays nearly constant for all values of As- As 
^5 1 the efficiency scales logarithmically in As- 

Figure [TOk shows the sliding power Vsuding/ls fo'^ the same base sliding values shown in Figs. [9] as functions of the 
slope angle, along with the rate of hydrodynamic work. For large sliding parameters As — > 1, the sliding power and 
rate of hydrodynamic work both decrease, but at the expense of a decrease in the swimming velocity. The decreasing 
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FIG. 9: The optimal infinite-length flagellar waveform with no bending cost {Ab ~ 0) but with non-zero cost of elastic sliding 
{As 7^ 0). (a) As the sliding cost becomes more dominant (increasing As) the slope of the resulting sawtoothed waveform 
decreases from 40.06° to a limiting value of approximately 33° for zero base sliding (ca =0). If the base sliding is larger 
the sliding cost is relatively less dependent upon the waveform, and the optimal solution does not change as significantly with 
increasing As- (b) The total and hydrodynamic efficiencies for ca = are labelled. The efficiency decays significantly with 
increasing As from the Lighthill value -q = 0.0858, but the hydrodynamic efficiency only decreases to a limit 0.0803. Total 
efficiencies corresponding to the ca constants in (a) are included as dotted lines. 



swimming velocity as a function of the slope angle is shown in Fig. [TOb . 
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FIG. 10: (a) The sliding power as a function of the slope angle |?/)| (degrees) for a selection of base sliding values ca, along with 
the rate of hydrodynamic work <I>. (b) The swimming velocity as a function of the slope angle, limiting to the Lighthill limit 
as As ^ (IVI « 40). 



IV. BODIES OF FINITE LENGTH 

The passage of periodic waveforms down along a flagellum of finite length introduces new degrees of freedom, 
namely vertical net motions, and rotations. The introduction of rotation can break time-reversal symmetry for bodies 
of non-half-integer numbers of wavelengths. In order to generalize the waveform for finite-length flagella, we include 
a wavelength parameter in the specification of the shape. Specifically, for finite-size swimmers, we optimize over the 
first n* Fourier modes describing the tangent angle: 

il){s) — a„ COS (2?™ k s) + 5„ sin {2Tm k s) , (47) 

n=l 
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where the wave-number k is to be determined as part of the optimization. Time is discretized into Tm uniformly 
distributed points on the domain t € [0, 1/k], and the body velocities xo(i) and 6(t) are determined at each time-step 
by inverting the three by three system, Eqns. [T] The number of time-steps must be sufficient to capture the activity 
of the highest Fourier modes in the travelling wave solution. We also insert an important constraint for optimal 
'swimming': we require 0{T) — 9{0) so that the organism does not rotate in circles over many periods. However, as 
we will show, this does not remove the possibility of a slow vertical drift perpendicular to the initial body orientation. 
We generally use here M = 2400 spatial gridpoints, Tm — 160 time-steps, and n* — 160 Fourier modes. The solutions 
reported here were checked against simulations using more refined spatial and temporal discretizations and greater 
numbers of Fourier modes, when possible. The results did not vary significantly with further resolution. 

For bodies of finite length there is an important degeneracy in the model in the pure hydrodynamic consideration 
of {Ab = As = 0). Given a body of finite length, the optimal solution must in fact be Lighthill's sawtooth function 
with infinitesimally small amplitude and infinitely many wavelengths. In this limit, there is no rotation, and hence 
there are no hydrodynamical costs associated with rotational work done on the fluid. Any non-zero bending costs will 
regularize the geometry, and a competition between body rotations and the bending costs associated with the number 
of wavelengths expressed by the body will ensue. In contrast, as we have shown above, the inclusion of a sliding cost 
(As > 0) does not regularize the optimal body shape. Hence, the optimal finite-length swimmer in the presence of 
a sliding cost must be the degenerate case of a body with infinitely-many wavelengths of infinitely small amplitude, 
with a slope as determined in the previous sections. We therefore consider below the influence of bending costs on the 
optimal finite-size flagellar waveform. The waveforms are not limited to any class of functions (other than periodic) 
and are determined by the numerical optimization. 



0.09 




total (rj, solid line) and hydrodynamic (ryn, dashed line) efficiencies. 

Figure [TT] shows the total and hydrodynamic efficiencies associated with the optimal shapes for Ab S [10~^, 1]. The 
optimal shapes for a selection of bending costs are also included, and are presented for a more direct comparison in 
Figs. [13 When the bending cost is very large compared to the hydrodynamic cost {Ab ~ 1), the optimal shape in an 
approximate sinusoid and expresses just beyond a single wavelength, k — 1.08. As the bending costs begin to decrease 
the optimal shape begins to express a slightly larger wave-number, k — 1.19 for Ab — 10^^. With further decreases in 
the bending costs the optimal shape takes on a sharper profile, and approaches a half-integer wave- number, k = 1.42 
for Ab = 10--*. 

While the shape appears to change continuously for bending parameters in the range Ab G [10^**, 1], we find a 
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FIG. 12: (color online) Optimal finite-size flagellar waveforms for a selection of bending costs Ab- As bending becomes less 
costly (decreasing Ab), the optimal shape expresses larger wave-numbers and sharper profiles, with a bias towards half-integer 
wavelengths. 



surprising transition between Ab — 10~^ and As = lO""*'^. While the efficiency appears to change continuously in this 
range of bending costs, the optimal shape jumps discontinuously to approximately the next half-integer wave-number. 
For Ab > 10^^ there is yet another discontinuous transition to the next half- integer wave-number, and the optimal 
shape becomes more and more like the infinite sawtoothed solution (77 = 0.858). We emphasize that this remarkable 
bias towards half-integer wave-numbers, which we observe over three full jumps, is an output of the optimization and 
not an assumed constraint. 

We have observed that while the efficiency appears to be continuous through the jumps in wave-number, it is not 
smooth (not shown here). The total efficiency decreases to a limiting value of 77 = 7.68 • 10^'* as the bending becomes 
exceedingly expensive {Ab 1). The hydrodynamic efficiency, on the other hand, decreases monotonically to a 
limiting value of rjH = 0.0606 for Ab 1. The hydrodynamic efficiency decreases by 30% in this limiting case, a 
more dramatic change than for the analogous body of infinite-length. Here the extra degrees of freedom, namely 
rotations and vertical drift, are more dependent upon body shape, and lead to larger variations in the dynamical work 
done to the fluid. Importantly however, this hydrodynamic efficiency of « 6% is still significantly above the efficiency 
of « 1% typically displayed by biological cells which utihze planar waves 

Other properties of the optimal finite-length fiagellum are shown in Figs. 1131 Data corresponding to the optimal 
shape are shown as solid points, but we also include hollow points to indicate values for certain locally optimal 
solutions. Figure [T3k shows the maximum curvature of the optimal shapes. Through the jump transition in wave- 
number there is a small jump in the maximum curvature, but the overall trend is preserved. As in the infinite-length 

— 1/2 

case, the maximum curvature of the finite-length body scales as approximately k ^ Ag a.s Ab ^ 0. However, 
particularly given the jumps in maximum curvature as the wave-number increases discontinuously the true asymptotic 
behavior as Ab ^ may not yet be well represented in this regime. As shown in Fig. 113b there is a distinct trend 
for decreasing bending costs towards half-integer wave-numbers, and there are also jumps to shapes of larger half- 
integer multiples at critical bending parameters. The left-right symmetry in shapes of half-integer wavelength serve 
to significantly decrease the body rotations throughout the motion. With smaller rotations, the body undulations can 
contribute more directly to lateral locomotion without performing much rotational work on the surrounding fluid. 

The swimming speed is shown in Fig. 113b . The locally optimal solutions with fc w 1.5 (shown in circles) give way 
to the globally optimal solutions with k w 2.5 (squares) at approximately Ab — 10^'' '^. The rate of increase in 
swimming speed for decreasing Ab becomes more rapid with the larger wave-number solutions, so that the swimming 
speed appears continuous but non-smooth. The swimming speeds are for all Ab smaller than the swimming speeds 
determined in the infinite-length consideration, as expected, since there are body rotations in the finite case which 
generally act to impede the lateral swimming motion. Figures [T3H and e show the maximum rotation angle ||0(i)||oo 
and the vertical shift (or slope of the velocity vector) as functions of the bending cost. Given sinusoidal or sawtooth 
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FIG. 13: (color online) Properties of the optimal finite-length flagellum as a function of the bending cost. The behavior 
limits to that of the degenerate, infinite-length solution as bending costs vanish {Ab ^0). (a) Maximum fiagellum curvature, 
(b) Wave-number, (c) Swimming speed, (d) Maximum body rotation, (e) Vertical body drift, (f) Wavelength-normalized 
amplitude. 



waveforms vi^ith integral numbers of wavelengths, it has been shown for small-amplitude motion that the rotation angle 
decreases with wave- number as ||0(i)||oo ~ 1/fc^ [30]. We have also found this scaling to hold for large amplitude 
waves (not shown). We observe non-monotonicity in the maximmn rotation angle for Ab ~ 10^^ and As ~ 10^^. 
The vertical drift also exhibits non-monotonicity in the same regions. The behavior near As = 1 is likely due to 
the transition from the extreme case of a single sinusoidal wavelength to the nearby (hydrodynamically preferred) 
half-integer wavelength. The small vertical drift of the swimming motion (also noted in (32|) is an effect which is third 
order in the wave amplitude for small amplitude waves, while the swimming velocity is second order in the amplitude; 
hence in small amplitude studies this drift is generally not observed. The drift decreases as the body takes on greater 
wave-numbers and undergoes smaller rotations. 

The hydrodynamic benefits of half-integer spatial modes is illustrated in Fig. [TH For fc « 1 the body experiences 
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a large rotation through the periodic motion, while the k « 1.5 mode for Ab = 10 shows damped rotations, and a 
more effective motion towards the left. The vertical drift is visible in the first case. 
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FIG. 14: (color online) Transition in the optimal waveform for a finite-size swimmer, (a) With Ab = 1, the increased cost 
of bending leads to a smoother waveform, and in turn to significant rotation. T — .994 is the fundamental period of this 
first shape, (b) For Ab ~ 10~^, bending is not as energetically costly, and a higher spatial mode is observed to be optimal. 
This corresponds to a reduction in rotation, and thus a more efficient motion opposite the direction of the travelling wave. 
Nearly-half-integral wave numbers benefit from their approximate left-right symmetry, which significantly decreases rotations. 
T = .994, as in (a), for comparison. 

Finally, Fig. [T3f shows the wavelength-normalized amplitude of the optimal finite-length flagellum, bk = 27r||Z||oo k. 
Perhaps surprisingly, even with the large rotations seen at small wave-numbers the optimal amplitude behaves very 
much like in the infinite-length case. For Ab = 1 is very large the wave amplitude is bk ^ 0.89, just larger than the 
infinite-length result. As the bending costs decrease the optimal waveform approaches approximately the same limiting 
amplitude seen in Fig. [7^, and for Ab — 10^'' we find bk ~ 1.315. Each jump to larger wave-numbers is accompanied 
by a jump in ||Z||oo. Hence, the optimal waveform appears to degenerate towards Lighthill's infinite-length sawtooth 
solution in a self-similar fashion. 

V. DISCUSSION 

In this paper, we have offered a physically-motivated derivation of the optimal flagellum shape. We have considered 
the optimal shapes of periodic, planar flagellar waves of both infinite and finite length in a model which, in addition to 
hydrodynamic dissipation, incorporates energetic costs of internal bending, sliding, and fiuid dissipation. For bodies 
of infinite-length, we have shown that the inclusion of a bending cost (or dissipation due to the presence of an internal 
fluid) regularizes the classical Lighthill sawtooth solution, and that the optimal waveform becomes very nearly (but 
not quite) a sinusoid. The inclusion of a sliding cost has been shown to decrease the amplitude of the optimal 
waveform, but the optimal shape is still a sawtooth with a jump in the slope at a finite number of points. For bodies 
of finite-length, we have shown that a degenerate solution, in which the body takes on infinitely many small amplitude 
waves, is regularized by the addition of any bending cost (or internal fiuid dissipation cost). With the exception of 
the case in which the bending is exceedingly expensive, the optimal shape has been shown to express approximately 
half-integer number of wavelengths, with a shape tending in a self-similar manner towards that of the infinite-length 
sawtoothed shape. This surprising result underlines the importance of minimizing the rotational work done on the 
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surrounding fluid during lateral swimming. In addition, for both the infinite- and finite-length cases, we have shown 
that the change in the hydrodynamic efficiency is relatively small, and remains well above the hydrodynamic efficiency 
of typical biological cells. 

The model presented here uses some simplifying assumptions, and leaves a number of open questions. First, the 
hydrodynamic description could be improved upon by the inclusion of non-local effects, for example using slender 
body theory, or a more complete three-dimensional method for thicker organisms such as nematodes. The work of 
Tarn jssi appears to indicate that the sawtooth form may not be regularized by the non-local fiuid interactions in 
the limit of zero bending costs, but that the number of expressed wavelengths may be decreased. Another exclusion 
in the work presented here is the possible presence of a head. Although most sperm cells have relatively small cell 
bodies (such as human spermatozoa), they can be large for some microorganisms and generally act to damp rotations 
imposed by the flagellar beating. In addition, the expression we used for the bending energy becomes invalid when 
the radius of curvature approaches the body radius. The formation of material or structural singularities has been 
considered by other authors, and this can also provide a barrier to the degeneracy mentioned above 




FIG. 15: Marine invertebrates spermatozoa, (a): Superimposed images of the headless spermatozoon of Lytechinus. (b): 
Spermatozoon of Chaetopterus exhibits non-integral spatial wave-numbers. Reproduced with permission by the Journal of 
Experimental Biology, Ref. p^ . 



In our opinion, the two most important implications of our study for the biophysics of swimming cells are the 
following. First, we have shown that a physically-motivated measure of internal elastic cost for the dcformable 
flagellum regularizes the hydrodynamically-optimal solution of Lighthill, and that this is done with only a small loss 
in hydrodynamic efficiency. Second, our results show the emergence of small numbers of wavelengths in the optimal 
solution when bending is at all costly (see Fig. Ilip . This result, which is likely to remain robust to improvements 
on the modeling such as the one discussed above, is consistent with observations on the morphology of eukaryotic 
flagella, in particular for spermatozoa (see the review in Ref. 21|). As an example, we reproduce in Fig. [15] the shapes 
of two marine invertebrates spermatozoa {Lytechinus and Chaetopterus) from Ref. [3]. In both cases, although the 
shapes are different from out optimal solutions, the presence of the half-integer wave-number morphology (fc « 1.5) is 
apparent. Our work constitutes therefore an attempt at a physical rationalization of this observed feature of eukaryotic 
flagella. We also observe that our optimal solutions display hydrodynamic efficiencies which are significantly above 
those of biological swimming cells, which are typically in the 1% range. Our solution could therefore also be considered 
as an appropriate 'initial condition' for further (more directly biological) optimization, at the expense of hydrodynamic 
efficiency. Finally, we note that another (less common) means of eukaryotic propulsion involves the passage of periodic 
helical waves down along the length of a flagellum. In this case, in addition to the costs discussed here, there may 
also be costs due to a twisting of the material, and dynein motors have been observed in some cases to exert twisting 
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moments on the axoneme 
considered in a future work. 



43|. The optimal shape of a helical flagellum under these energetic constraints will be 
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